

Hi_UCI<-c(2.0858255, 2.050654,    2.0154918,   1.9803396,  1.9451978, 1.910067,  1.8749477,  1.8398407,  1.8047465, 1.7696659 ,1.7345997,  1.6995487, 1.6645139, 1.6294962, 1.5944967, 1.5595165, 1.524557,  1.4896194, 1.4547053,  1.4198162,  1.3849539, 1.3501203,  1.3153175,   1.2805479,  1.2458138, 1.2111182 )

Lo_UCI<-c( -.00518132,  .01426948 ,   .03317144 ,  .05157636 , .06953352 , .08708885 ,  .10428447, .12115858, .13774552, .15407594, .17017707, .18607304 , .20178519, .21733235, .23273118, .24799637, 

					.26314091, .27817632, .29311279, .30795937,  .32272411 ,  .33741419, .35203602, .36659532 , .38109724 , .39554639                   )


effect<-.79

deviations<-seq(.26,.01,by=-.01)

prop_dev<-deviations/effect


effect<-rep(effect,length(prop_dev))



pdf(file="Figure_A12_left.pdf")

plot(prop_dev,effect,pch=16,cex=.5,ylim=c(-.02, 2.1),xlab="Prior on %  Direct Effect Exclusion Violation",ylab="Effect",main="Within")

for(i in 26:1){

lines(

c(prop_dev[i] , prop_dev[i]  ) , c(Lo_UCI[i], Hi_UCI[i]) ) 



}


abline(h=0,lty=2)

dev.off()